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Abstract 



In recent years, a large number of papers have appeared that dealt with e + e~ 
pair production in heavy ion collisions at high energies. To a large extent these 
studies were motivated by the existence of relativistic heavy ion accelerators all 
over the world. There pair production can be studied in so called "ultra-peripheral 
collisions" , where the ions do not come close enough to interact strongly with each 
other. Various different methods have been used and it is the purpose of this review 
to present a unified picture of the present status of the field. The lowest order Born 
result has been known for more than seven decades. The interest and focus is now on 
higher order effects for values of Za ^$ 1, where Z is the charge number of the ion. A 
similar problem appears for the Bethe-Heitler process, the production of e + e~ pairs 
in photon-nucleus collisions. It was solved essentially some five decades ago by Bethe 
and Maximon. The result of Bethe and Maximon can also be recovered by summing 
over a class of Feynman diagrams to infinite order. These results can be used for a 
study of Coulomb corrections in nucleus-nucleus collisions. Indeed, the major part 
of these corrections have a structure closely related to the Bethe-Maximon solution. 
There are additional terms which give a small contribution to the total cross section 
at high energies. Their importance can be enhanced by concentrating on small im- 
pact parameters. An interesting exact solution of the one-particle Dirac equation in 
the high-energy limit was found independently by several authors. This led to some 
discussion about the interpretation of these results within QED and the correct 
regularization necessary to get the correct result. The dust of previous debates has 
settled and, indeed, a consistent picture has emerged. Another interesting higher 
order effect is multiple pair production, which we also discuss. We compare exper- 
imental results obtained recently at RHIC for free and bound-free pair production 
with theoretical results. We also make some more remarks on the physics of strong 
electric fields of longer duration. A new field is opened up by ultra-intense laser 
pulses. We argue that due to the short interaction time in ultraperipheral heavy ion 
pair production can be well understood in the frame of QED perturbation theory. 

Key words: Electron-Positron Pair Production, QED of strong fields, Coulomb 
Corrections, Multiple Pair Production, Relativistic Heavy Ion Collisions, 
Ultraperipheral Collisions. 



Contents 

1 Introduction, overview and purpose 

2 Lepton pair production by a high energy photon in the strong 



Coulomb field of a heavy nucleus 



11 



3 



Lepton pair production in relativistic heavy-ion collisions using 
Feynman diagrams 







2 



4 Exact solution of the one-particle Dirac equation in the 
ultrarelativistic limit and what one can conclude from that 

5 Coulomb corrections for the heavy ion case, based on the 
Bethe-Maximon approach [22] 

5.1 Coulomb corrections in the equivalent photon approximation: the 
approach of Ivanov, Serbo and Schiller [2^ 

5.2 How to regularize the exact solution of Baltz and McLerran and 
Segev and Wells: the approach of Lee and Milstein \2A 

6 Higher Order Effects in Electron Pair Production: Multiple Pair 
Production and Bound-Free Pair Production [2I 

6.1 Multiple pair production 

6.2 Bound-free pair production |34 

6.3 Higher order corrections to muon pair production |35 

7 Transition from the adiabatic to the sudden regime [37 

8 Comparison to experiment and outlook to LHC 

8.1 e + e~ pair production in ultraperipheral collisions at RHIC [4(] 

8.2 Bound-free pair production at RHIC and the forthcoming 
LHC(Pb-Pb) [H 

8.3 Outlook on LHC E 

9 Conclusion |45 
References 47 



3 



1 Introduction, overview and purpose 

Quantum electrodynamics (QED) is the best-tested theory we have [TlTSllBTf] . 
It is a quantum field theory which describes the gauge invariant interaction of 
charged particles with photons. With this theory one is able to describe high 
precision experiments like the Lamb shift and the magnetic moments of the 
electron and muon with great accuracy. In the present review we are dealing 
with the fundamental process of particle production in high energy scattering 
processes. Experiments up to now have been not very accurate in their results. 
Also, the theoretical precision was not very high. This is not due to principal 
problems of the underlying theory, the QED. It is more related to technical 
difficulties, as we shall see in the following. An especially interesting process 
is the production of electron-positron pairs under various conditions. 

Positrons were first observed in cosmic ray interactions in 1932 [5]. The the- 
oretical calculation of e + e~ pair production in fast nucleus-nucleus collisions 
goes back to the early days of QED [Bj. The full result in lowest order pertur- 
bation theory (Born approximation) is due to Racah [7] . These results are also 
discussed in [8]. In those calculations the pair production is due to the colli- 
sion of a pair of (virtual) photons. The theoretical, as well as, experimental 
interest in these processes was spurred in the last decades by the relativistic 
heavy ion facilities AGS at Brookhaven, SPS at CERN and more recently 
RHIC at Brookhaven and LHC at CERN. It was realized that the old lowest 
order answers are no longer sufficient for the strong Coulomb field surrounding 
the heavy ions at these high energies. There are (new) sizeable and interesting 
higher order effects. 

The subject of this review is precisely the situation envisaged in these early 
works of the thirties of the past century, i.e. pair production in high en- 
ergy nucleus-nucleus collisions. The theoretical treatment now concentrates on 
higher order effects. One type of effects is the influence of higher orders on the 
one-pair production. Very similar "Coulomb corrections" have been discussed 
already before in pair production by a single (real) photon in the Coulomb 
field of a nucleus. As we will see a lot of the experience gained there can also 
be applied to the case of two heavy ions. Another — new and more spectac- 
ular — effect is the production of many pairs in a single collisions, clearly a 
higher order effect. Various groups have been working on these problems in 
the past two decades using various methods and approximations. The main 
purpose of this review is to discuss these various methods and to establish 
the relationship of the different approaches. Since QED is a well established 
theory, such a connection is to be expected. 

Before dealing specifically with our subject of relativistic heavy ion collisions, 
it is appropriate to put the present subject of e + e~ production in relativistic 



4 



nucleus nucleus collisions into a wider perspective: slow collisions of heavy 
ions of charge numbers Z\ and Z 2 are a way to create, at least for a short 
time, atoms with a charge United atom = Z x + Z 2 - It is of special interest to 
study collisions where the sum of the charges United atom > 137 (173). In 
these collisions so called "overcritical fields" are reached for some time where 
spontaneous e + e~ production is possible, see [9] or the detailed review of this 
field in [ID] . The supercritical field is reached for a charge Z = 137 in the 
case of a point nucleus, where one may say that the lsi/2-level dives into the 
negative energy continuum. For extended nuclei with a finite radius, this limit 
is reached at about Z = 173 [DJ. 

The paradigm of e + e~ pair production is the Schwinger mechanism: the pro- 
duction of pairs in a constant spatially uniform electric field [11] , see Fig. [IJ To 
reach a sizeable rate for spontaneous production of e + e~ pairs a field strength 
of the order of or above a certain critical field strength E c is required. The 
rate of (one) pair production per unit volume and time for a given constant 
field E is given by 



d 4 n e+e - c . E c 

jexp(-7v— ), 



d 3 xdt 4vr 3 A 4 v E 



where the critical field strength is given by 

2 S 

E c = — — ~ 1.3 x 10 16 Wcm. 

en 



The electron mass is denoted by m and its Compton wavelength is given 
by X = — . This formula clearly shows the non-perturbative nature of the 

J mc J 1 

problem: there is no power series expansion of the e + e~ production rate in 
terms of the electric field strength E. With new developments in laser physics 
it seems that this field strength can be reached, e.g., at the FEL at DESY and 
the interest has grown enormously, both in theory and experiment [I5fl%] . In 
contrast to relativistic heavy ion collisions, the created fields in this case are 
of a comparatively long duration. This is made more quantitative in [T5] and 
will be dealt with in Sec. [7] 

The main purpose of relativistic heavy ion collisions is the production and 
observation of the quark gluon plasma which yields information about QCD 
at high temperature and densities. It has gradually been realized over the last 
decades that the very strong electromagnetic fields present in distant heavy 
ion collisions are also useful: they provide a strong flux of equivalent photons 
for photon-hadron and photon-photon physics with hitherto unknown ener- 
gies. The term "ultraperipheral collisions" (UPC) has been coined to describe 
collisions where the impact parameter is larger than the sum of the two nu- 
clei: in this case only electromagnetic interactions occur between the ions. 
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V(x) 




Fig. 1. Pair production in a strong static field E. The total potential is given by the 
potential energy (solid line) of the electron in the binding potential V{x) (dashed 
line) and in the static electric potential eE. We may consider the electron bound by 
a potential of depth ~ 2mc 2 . In the Schwinger process the electron tunnels through 
the barrier and an e + e~ pair is created, see [12]. 



Fig. 2. The electromagnetic fields of two heavy ions in an ultraperipheral col- 
lision. These fields can be decomposed into a spectrum of quasireal photons 
(Weizsacker- Williams approximation). The collision of two photons gives rise to 
a dilepton pair. The collision of photons radiated from each nucleus are a copious 
source of e + e~ pairs. 

One special case of photon-photon physics in UPCs is e + e~ production, see 
Fig. [2j Due to the smallness of the electron mass e + e~ pair production is very 
copious, it is of the order of kilobarns under RHIC and LHC conditions. The 
physics of "ultraperiphal collisions" is reviewed in [T6]I17M18|[T9] and [20]. A 
recent workshop at ECT*/Trento, Italy was devoted to the subject of UPC 
[2~T] . The status of e + e~ pair production was also discussed there, see [22] . 

Let us make some order of magnitude estimates of the fields acting in the case 
of relativistic heavy ion collisions: the electromagnetic fields encountered in 
these collisions are extremely strong, they are of the order of 



where we assume that the ions move on straight line trajectories with impact 
parameter b and 7 is the Lorentz factor in the collider system. With the values 
of Z = 79, b = b min = 15 fm and 7 = 100, which are appropriate for RHIC, we 




E, 



Ze^y 



(3) 



max — 
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find a maximum field strength of E max = 4.9 x 10 16 V/cm. For LHC conditions 
(Z = 82,b min = 15/m,7 = 3000) we get £ max = 1.5 x 10 18 V/cm. These are 
the values viewed in the collider frame. Viewed from the other nucleus the 
corresponding Lorentz factor is 7 ion = 27 2 — 1 and the electric fields are even 
larger by a factor of 7. 

These field strengths are comparable to the critical field E c for the Schwinger 
mechanism. However, these fields act only for a very short time 

At ~ — , (4) 

7f 



where v ~ c is the ion velocity. For b = X ~ 386fm one obtains the values of 
At ~ 10~ 23 s for RHIC and At ~ 3 x 10 _25 s for LHC conditions, respectively. 
The momentum transfer Ap given to a test particle with a charge e along a 
path with impact parameter b is independent of 7. For v ~ c and b = X it is 
only modest 

Ze 2 

Ap ~ eE max At ~ — ~ Zamc. (5) 



This small value gives a strong hint that the use of perturbation theory is still 
appropriate in contrast to the situation of the Schwinger mechanism mentioned 
above. Surely there are higher order effects, but they can be dealt within 
higher order perturbation theory, see, e.g., [23]. One goal of this review is 
the explanation of the results for the effects beyond lowest order perturbation 
theory. An obvious question is about the transition from the slow to the fast 
collisions. 

In section [2] we recall a simpler but intimately related problem: the Bethe- 
Heitler process, i.e. the photoproduction of e + e~ pairs in the nuclear Coulomb 
field. It shows already many of the characteristic aspects of the heavy ion case. 
Higher order effects were treated in [24T25] by using the Sommerfeld-Maue 
wave functions: these wave functions are appropriate high energy solutions 
of the Dirac equation and thus take higher order effects into account to all 
orders. This is a textbook example (see e.g. [26]). It was shown in [27] how 
this approach is related to the usual Feynman graph technique. 

In section [3] we come back to the case of ultraperipheral heavy ion collisions. 
We first review the Feynman graph approach to pair production in relativistic 
heavy ion collisions. The lowest order (Born approximation) cross section was 
given by Landau and Lifshitz [6] and Racah [7] in the thirties of the past cen- 
tury. In contrast to the photoproduction on a nucleus one has a well-defined 
impact parameter in heavy ion collisions and impact-parameter dependent 
probabilities have been studied as well. Higher order effects were studied in the 
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Fig. 3. e + e~ pair production is classified by the number of photons attached to each 
ion. We distinguish for classes (i) n = n' = 1 (ii) n = 1, n' > 1 (iii) n > 1, n' = 1 
and (iv) n > 1, n' > 1. 

last decades, mainly by the russian schools. We note here that due to the high 
mass of the projectile bremsstrahlung pair production is negligible at the high 
energies, as explained, e.g., in [26], see also [2H]. Bremsstrahlung pair produc- 
tion is relevant for lower beam energies, where the (timelike) bremsstrahlungs- 
photon is converted to an e + e~ pair. 

As we will see in the following we have an important classification of the terms: 
We denote the number of photons attached to nucleus 1 by n, to nucleus 2 by 
n and distinguish four classes, see Fig. [3j 

(i) n = 1 and n' = 1 

This is the lowest order treated by Landau and Lifshitz [6] and Racah 

[7]. The cross section increases with the Lorentz factor 7 like (log7) 3 and is 
the dominant term at high energies. 

(ii) n = 1 and n' > 1 

(iii) n > 1 and n! = 1 

These terms are the main correction and well understood as we will show. The 
problem is essentially the same as the one for the Bethe-Heitler process, and 
it is solved in an analogous way. 

The case, which is specific to heavy ion collisions and which contains all the 
problems is the last one: 

(iv) n > 1 and n' > 1 

This case cannot be considered to be solved, but at least it is reassuring that 
one can show that it is only a small contribution to the cross section. This 
classification is directly evident in the Feynman diagram approach of Sec. [31 
but it also plays the major role in the approaches of Sec. [51 When the ion 
velocity decreases or for small impact parameter the Feynman diagrams of 
class (iv) will become more important. So collisions at intermediate and small 



energies become increasingly difficult to treat theoretically. Molecular orbital 
approaches will then become relevant. This is quite analogous to the situation 
in atomic physics for the excitation of atoms by ions in ion-atom collisions. 
This is an interesting topic but it is outside the scope of this review. 

In Sec. H] an analytical solution to the (one-particle) time dependent two- 
center Dirac equation in the high energy limit is discussed. The authors of 
[29|30|31|l32] were led to the (erroneous) conclusion that the cross section cor- 
responding to the exact solution is the same as the cross section in the Born 
approximation. But this is in obvious contradiction to the special case of the 
Bethe-Maximon corrections to the Bethe-Heitler formula. There are two as- 
pects in this work: one is related to the interpretation of the one-particle Dirac 
equation to the corresponding (many-body) Feynman diagrams describing pair 
production. We will show how the two can be related, see, e.g. [3U1IM] . The 
other aspect has to do with the proper regularization of the expression, which 
was used by [29.31,32.35] in the high energy (sudden) limit. 

In Sec. [5] we discuss the work of Ivanov, Schiller and Serbo [36] and Lee and 
Milstein [37J. The latter work is based on the analytical solution of the Dirac 
equation in Sec. HI In Sec. E] we discuss some aspects of class (iv): two- or 
more photons emerge from either nucleus: in addition to the 'Coulomb correc- 
tions', which are small, there is a new genuine higher order effect: the multiple 
pair production in a single collision. Another higher order effect (belonging 
essentially to classes (ii) and (hi)) is bound- free pair production: the electro- 
magnetic interaction of the ion (to all orders) with the electron leads to a 
bound state (i^-,L-shell,. .. capture). This is a small but interesting fraction 
of the total pair production, the cross section scales only with log(7), instead 
of (log(7)) 3 for free pairs. In Sec. [7] we make some remarks about the transi- 
tion from the non-perturbative slow to the perturbative fast collisions. This 
is a difficult question. We discuss especially the results of [T5]. An interpola- 
tion between the two limits is achieved analytically for the case of spatially 
constant electric fields. In Sec. |8] we give a comparison of theory to experi- 
ment: in the last five years since our previous review [18] RHIC has come into 
operation. We compare the experimental results of STAR on e + e" pair pro- 
duction in UPC [38J to lowest order QED calculations and discuss the results 
of bound-free pair production at RHIC. Accelerator physicists are anxious to 
test the theory of bound-free pair production because this process is crucial 
for the operation of the forthcoming LHC(Pb-Pb). The bound-free capture 
process will be the ultimate limit for the Pb-Pb luminosity. Conclusions and 
an outlook are given in Sec. EJ 

Throughout this review we use H = c = 1 and a = e 2 , unless it seems of some 
pedagogical value to include these constants explicitly. We also assume rela- 
tivistic high energy collisions, which are characterized by 7 3> 1. We assume 
that even log(7) 3> 1. This is well fulfilled for RHIC and LHC conditions, less 
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for AGS and SPS energies. 



2 Lepton pair production by a high energy photon in the strong 
Coulomb field of a heavy nucleus 



In this chapter we recall the well known theoretical treatment of the Bethe- 
Heitler process, the e + e~pair production by a photon in the Coulomb field of 
a nucleus. This is simpler than the problem of pair production in the field of 
two nuclei, but shows already many of the important aspects of the nucleus- 
nucleus case. It corresponds to the cases (i), (ii) and (iii) (see Sec. [TJ, where 
only one photon is attached to one of the nuclei. Case (iv) is the one, which 
is specific to the nucleus-nucleus problem. 

The Born approximation result (case (i)) is derived in textbooks using the well 
established tools of QED, see e.g. [26]. The celebrated Bethe-Heitler result for 
the total cross section of the process 



where the photon energy is denoted by u. The Coulomb field of the nucleus is 
treated as an external field. It is a characteristic feature that the cross section 
rises logarithmically with the photon energy. Screening due to the atomic 
electrons leads to a constant cross section for high energies. This screening 
effect is of little interest for us here, as we are mainly interested in the case of 
bare ions at the colliders. 

The higher order effects in pair production at high photon energies and also 
for large energies of electron and positron compared to their rest mass, were 
studied by Bethe, Maximon and Davies [2^1123] . 

Their results are also described in detail in many textbooks (see especially 
|26j). For the derivation of the cross section these authors used Sommmerfeld- 
Maue wave functions for the electron and the positron, respectively. These 
wave functions are approximate solutions of the Dirac equation in the nuclear 
Coulomb field valid for high energies and small scattering angles. Thus they 
include, in a certain way, higher order effects up to all orders. They are given 
for the electron with energy e_ and momentum by: 



1 + Z^e + + e- + Z 



(6) 



in the case of no screening is given in the high energy limit by 




(7) 




=e ip -" r 




ct • V 



)F(—iZa, 1, — i(p_r + J?_ ■ ~r y ))u(p-)(8) 



v 7 ^ 



2e_ 
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with the normalization 

iV_ = exp(^Zct)r(l + iZa) 



(9) 



depending on the T-function and the confluent hypergeometric function F. 
The Dirac spinor for the electron in the free case is given by u and ~ct are the 
Dirac matrices. 

Similarly the wave function for the positron is given by: 

^+J +p+ = -^ e -^+' y (1 + i—^-)F(-iZa, l,i(p + r + f + • T))u(p+)(10) 



with 



N + = exp(--Za)r(l - iZa). (11) 



Inserting these wave functions into the standard matrix element, differential 
and total cross section can be calculated analytically, leading (after neglecting 
some small terms in the high energy limit) to the expression for the total cross 
section: 

28 a 3 Z 2 „ 2uj 109 „.„ 
9 vnr m 42 



where the Bethe-Heitler cross section, see Eq. (|7J above is now modified by a 
term containing the function f(aZ). It is defined by: 

OO 1 

f(Za) = lE + ite¥(l + iZa) = (Za) 2 £ (13) 

™(™ 2 + (^«) 2 ) 



where 7^ = 0.57721... is the Euler-Mascheroni constant and \P the loga- 
rithmic derivative of the T-function. Note that in this way higher order effects 
are established without making use of perturbative QED. In fact, Feynman 
graphs up to infinite order are effectively summed up by the use of these 
Sommerfeld-Maue wave functions. 

It is very interesting to see that the same results can also be derived directly 
(among other things) by the standard method of summing up all Feynman 
graphs [27], see Fig. [U This establishes the equivalence between the standard 
higher order QED approach and the calculations using Sommerfeld-Maue wave 
functions. This is very important for our problem below: As discussed in the 
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(b) 



Fig. 4. (a) Pair production in photon-nucleus scattering in lowest order of Za. 
(b) Higher order term. These terms are summed up in the Bethe-Maximon theory 
using the Sommerfeld-Maue wave functions, or by summing the Feynman graphs 
according to the approach of [2.7]. 

following chapters, there are again the Feynman diagram approaches to the 
nucleus- nucleus problem, and the 'Sommerfeld-Maue type' of approaches. 

In the following, we very briefly describe the work of Ivanov and Melnikov |27j . 
A typical type of Feynman graphs to be calculated is shown in FigJH (which 
is Fig. 1 of }27j). For N — 1 this corresponds to the Born approximation: the 
e + e _ pair is produced by the collision of a real photon (which can also be 
virtual in [27]) and a virtual photon coming from the nuclear Coulomb field. 

Many steps are needed to sum up the Feynman graphs in the high energy limit. 
We sketch some of the important points: The basic kinematics is pi + P2 = 
Pi + Qi + <?2 and it is convenient to regard the nucleus as a light particle with 
mass m and charge Z. We assume 2pi -p 2 = s^> m 2 , Q 2 . We denote p\ = —Q 2 ; 
for a real photon we have Q 2 = 0. There are two important scales for the 
square of the transverse part of the transferred momentum q = q\ + q 2 — p\\ 
q 2 ~ m 2 , Q 2 and q 2 ~ m 4 / s, Q 4 / s. The second low momentum transfer region 
is only important for the case where one photon is exchanged, it is responsible 
for the logarithmic rise of the cross section. For the case of more than one 
photon exchange with the nucleus the first region is the dominant one. 

For the high energy limit considered in [27] the so-called Sudakov variables are 
used as the appropriate variables. They are light-like vector p~\ and p 2 , which 
are pointing almost in the direction of p\ and p 2 . They are given by 

Q 2 m 2 

Pl=Pl-\ P2, P2=P2 Pi- (14) 

s s 

These Sudakov variables can also be related to the so-called light-cone vari- 
ables. The use of these variables allows to express the four- dimensional loop 
integrals occuring in the calculation of the Feynman diagrams in terms of the 
transverse momentum only. The transverse momentum integrals themselves 
are expressed in coordinate (impact-parameter) space and the high energy am- 
plitude contains the so-called impact factors. These are scalar and vector struc- 
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ture quantities for the exchange of iV-photons. Using recursion relations and 
other tricks to which we can only refer to the original paper [27] the authors 
are able to sum them. 

The full result for the cross section for the pair production due to a virtual 
photon 7* and where all the graphs with N = 1 to N = oo are summed is 
given in Sec. IV B of [27J . For a real incident photon the result was found to be 
identical to the one obtained previously by using the Sommerfeld-Maue wave 
functions. The methods of [27j can also be applied to the relativistic heavy 
ion case, which is the subject of the next section. 

Just as in [27], we are especially interested in the high energy limit, but it is 
also of more general interest to study corrections in K For this question we 
refer to Ref. [4"0] . 

We note that the theoretical results in the high energy small scattering angle 
limit are symmetric with respect to the exchange of e + and e~. It is certainly 
not the case for bound-free pair production, a process to be discussed in Sec. El 
see [4~Tf4~2] . We also mention without further discussion that another asymme- 
try was introduced by Brodsky and Gillespie [43J in their second order Born 
treatment of large angle e + e~ pair production. 
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3 Lepton pair production in relativistic heavy-ion collisions using 
Feynman diagrams 



In this paragraph we want to briefly review the work on e + e~ pair production 
in relativistic heavy ion collisions based on the summation of Feynman graphs 
in the high energy limit. The problem is formulated in [Hj. The process to be 
studied is 

A( Pl ) + B(p 2 ) - e+(g+) + e"(g_) + A(p[) + B(p' 2 ), (15) 



where the corresponding four-momenta are given in the brackets. 

In the paper [33] the first terms in the amplitude of e + e~ pair production in 
the Coulomb field of two relativistic heavy ions are calculated. The Sudakov 
technique is used which simplifies the calculations at high energies, see Sec. |2j 
In lowest order in a (corresponding to class (i), see Fig [3]) the classic result of 
Racah is obtained 



a = ^2M^(L 3 _ 2.2L 2 + 3.84L - 1.636), (16) 
Tim 2 27 



where L = ln^^) = hi(7 2 ) for symmetric collisions. For details we refer to 
reference [33] ■ Also further references are given there, which actually give the 
steps needed to derive the Racah formula by their methods. 

Next they consider Coulomb corrections. The first correction corresponds to 
the case where one nucleus emits one and the other nucleus emits two photons 
(n = 1, n' — 2), corresponding to the lowest order diagram of class (ii). This 
contributes to the Coulomb corrections with the pair produced in an C-odd 
state. Altogether six Feynman diagrams contribute to this amplitude. 

The first C-even correction involves graphs where n — l,n' — 3 (belonging 
again to class (iii)) and n = 2, n! = 2 (the lowest order diagram contributing to 
class (iv)). Whereas the former case proved to be tractable analytically, obsta- 
cles were found for the second case. An explicit calculation of the n = 2,n' = 2 
case is reported in [33]. Their result suggests that the familiar eikonalization 
of Coulomb distortions breaks down for the oppositely moving centers. 

Lepton pair production in relativistic ion collisions to all orders in Za with 
logarithmic accuracy is studied in [46] . The amplitude is separated into the 
four terms 

M = M {i) + M {ii) + M {iii) + M {iv) (17) 
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corresponding to the four different classes defined in the introduction. Different 
terms in the absolute square of the amplitudes can then be classified according 
to their energy dependence, i.e. by their dependence on different powers of 
the large logarithm L = 111(7172) = ln(7 2 ). The absolute square of the Born 
amplitude |M(j)| 2 leads to the famous rise of the cross section with L 3 . It 
should be noted that the Racah formula (see Eq. (fl6l) ) contains also terms 
proportional to L 2 and L. As discussed in jH] the next powers of L come 
from the absolute square of Mui) and Mu^ and their interference with the 
Born amplitude May The leading term of their result coincides with the one 
found by other methods, which will be described below. Their result for the 
Coulomb corrections in order L 2 is of the 'Bethe-Maximon type': 

a c = 2 ^ Z2l f 2 {f{Z ia )f{Z 2 a)L 2 + O(L)), (18) 
where f(Za) was defined in Eq. ([131) . 

A direct manifestation of strong fields is multiple pair production. Early work 
on this subject started with the observation that the impact parameter depen- 
dent total pair production probability calculated in lowest order perturbation 
theory is larger than one, i.e. P^'(b) > 1. More on multiple pair production 
will be given in Sec. [6j The analysis of [17] fits in naturally within the approach 
to study Feynman graphs in the high energy limit. The general Feynman di- 
agram is shown in Fig. [51 It has N s light-by- light (LBL) blocks N e photon 
exchanges between the ions Z\ and Z 2 (not present in an external field ap- 
proach) and N pair production processes. Coulomb corrections are neglected. 
They assume Z\ )2 ^> 1, Z\a ~ Z 2 a < 1 and Z\Z 2 a > 1. The probability of N 
pair production is shown to obey a Poisson distribution. 

In [18] a new approach to lepton pair production in the Coulomb field of 
two relativistic nuclei was developed. It was shown that for certain sums of 
finite terms of the Watson series cancellations lead to infrared stability of 
the amplitude (i.e. the result does not depend on a screening radius of the 
Coulomb field). 
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Fig. 5. A general diagram for the process of iV-pair production in ultraperipheral 
nucleus-nucleus collisions is shown. It contains N e photon exchanges between the 
nuclei, N s light-by-light scattering blocks and N e + e~ pair production diagrams. 
The symbols for light-by-light scattering blocks and pair production diagrams are 
given in terms of the Feynman graphs in Fig. (b) . 
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4 Exact solution of the one-particle Dirac equation in the ultrarel- 
ativistic limit and what one can conclude from that 



An interesting topic is the solution of the (single-particle) Dirac equation 
for two countermoving ions and how this relates to the process of pair pro- 
duction. This was strongly triggered by the observation, that in the limit of 
ultrarelativistic nuclei and using an appropriate gauge, the expression for the 
electromagnetic interaction simplifies and the Dirac equation can be solved 
analytically in a closed form |49j . 

Before we discuss this specific solution of the Dirac equation, it is useful to 
remember that e + e~ pair production is a problem in relativistic quantum 
field theory. It is inherently a many-particle theory from the beginning, where 
particles can be created and destroyed. In contrast to this situation the familiar 
problem of ionization of an atom by a non-relativistic heavy projectile can 
indeed be dealt with by the time-dependent Schrodinger equation: the atomic 
electron (say in a hydrogen atom) is described by a Hamiltonian Hq = T+e 2 /r. 
The classical motion of the projectile causes a time-dependent perturbation 
V(t) and the problem is to solve the time-dependent Schrodinger equation for 
the wave function $ of the electron 

z^r = (H + V(t)Mt). (19) 



This equation can be solved by a number of well-established methods like 
perturbation theory, see e.g. [50]. For special cases, like the harmonic oscillator, 
there are also full analytical solutions. Alternatively, one can solve numerically 
the Schrodinger equation for the wave function $(£) and project on specific 
asymptotic states in order to determine the cross-section for transitions to 
these final states. 

In the case of e + e~ pair production Eq. (fT9|) has a different meaning: $ is now 
a state vector, the term V(t) = j^A^ is the interaction of the time-dependent 
(external) field A^{t) with the electromagnetic current operator = ^7^. 
This current can be written in terms of electron and positron creation and an- 
nihilation operators. One gets different terms describing electron and positron 
scattering, as well as, e + e~ pair production and annihilation. This problem is 
solved by the Dyson expansion (see, e.g., [26] or other textbooks on Quantum 
field theory). This approach then leads to amplitudes for e + e~ pair production 
not only for single pair production, but for N = 0, 1, 2, . . . pair production. 
Multiple pair production will be discussed later in Sec. El Here we want to 
discuss only the relation between a solution of the Dirac equation and the 
pair production process. That there is a problem with the interpretation of 
the one-body Dirac equation can already be seen in the discussion of the free 
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Dirac propagator, see e.g. Chap. 2-5 of [12]. The solution of the (free) Dirac 
equation (with V(t) = 0) can be written as 

*(i 2 ,f 2 ) = y"rf 3 x 1 fs:(t2,x 2 ;ti,fi)7o^(ti,fi), (20) 



where K(x2,Xi) is the retarded kernel. However, as is remarked on p. 90 of 
[12] , the hole theory suggests the introduction of a different propagator, the 
Feynman propagator. It appears in a natural way in the quantized field theory. 

In our case here, that is, for the interaction with only an external field we can 
find at least a relation between the two pictures, which we quickly want to 
sketch here. For further details we refer to [3"3"llo"T] . Only for certain processes 
one can neglect the many-particle aspect. This is for example the case for 
photoproduction, where the interaction with only one photon and a static 
Coulomb field leads to only one electron lifted from the negative sea to the 
positive continuum. 

Baltz, Gelis, McLerran and Peshier discuss the meaning of the exact solution of 
the Dirac equation in Ref. [30] . They make a valid point: the interpretation of 
the solution of the one-body Dirac equation is not straightforward. They show 
that the total cross-section obtained in this approach does not correspond the 
exclusive cross-section to produce exactly one pair, but to one where multiple 
pair production cross section are also included (each weighted by the number 
of particles produced). This is rederived in a very straightforward manner in 
[33]. 

The relation between the two approaches can be found in the following way, 
see [23] or [2152"] for details. The 5-matrix of the Dirac equation relates the 
annihilation operators for the electrons b and creation operators for positrons 
d, respectively, in the initial (i) and final (/) state in the following way: 




Here the indices + and — refer to the positive or negative energy, that is S-+ 
is the transition from a positive to a negative energy state. It is this matrix S 
that is calculated by solving the Dirac equation. 

For the Feynman boundary condition one prescribes the initial electron and 
the final positron (following the picture of Feynman that the positrons are 
electrons traveling back in time) and solves the field equations to get the final 
electron and initial positron states. For single pair production we have one 
positron in the final state and calculate those processes leading to an electron 
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in the final state. Therefore we need to rewrite the relation in Eq. (f2~Tj) as one 
between the b % and d-f + on the one and the and d %+ on the other side. One 
gets 

bf = (S ++ - + S + ^Szid f+ . (22) 



The first term is of interest for electron scattering, whereas the second term is 
the one we are interested in. In addition we need to add the vacuum-vacuum 
transition amplitude C, see Sec. below, to get for the matrix element for the 
single pair production process 

M kl = (0\b{d{\0) = C(S + -S^_) M . (23) 



Having a complete solution of the Dirac equation one can now in principle 

invert S (the amplitude for the scattering of positrons to positrons) in order 

to get the correct single-pair production amplitude. 

In the papers of Segev and Wells [3~Tf3"2] . Baltz and McLerran [29] and Eich- 
mann, Reinhardt, Schramm and Greiner [35] a remarkable analytic solution of 
the (one-body) Dirac equation is found. It was shown subsequently by Lee and 
Milstein how to regularize this expression and to find a useful solution for the 
Coulomb correction problem (see the next Sec. [5]). By a suitable gauge trans- 
formation [29] the potential created by the counter-moving ultrarelativistic 
ions can be written as 

V(p, z, t) = -6(t - z)a ■ VA~(p) - 6(t + z)a ■ VA + (p), (24) 



where 



A ± W = -Z„1»(*±W!. (25) 



Except at x ± = -j^(t±z) the electron or positron propagates as a free particle. 
As is remarked in [29], to construct the propagator for the electron, one needs 
to solve a boundary value problem with free propagation everywhere except 
at the surfaces of discontinuity at the light cone x = 0, see Fig. [6] The 
final result is their Eq. (51). We do not give it here but refer to the next 
section, where the work in [53] and [51] is discussed. There one starts from 
their formula Eq. fl24l) . which is now contained in Eqs. fl30|) - fl33l) of Sec. [5] and 
show how to regularize this result correctly. 
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/ (z+t) 



E<0 
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Fig. 6. The special structure of the electromagnetic interaction can be seen in this t-z 
plot. Due to the Lorentz contraction, the electromagnetic fields are localized in two 
sheets corresponding to 2 = ±t. In the 'retarded' or 'Dirac sea' approach (left), an 
electron with negative energy comes from t = — 00, crosses the field of each ion only 
once before leaving as an electron with positive energy. In the 'Feynman' approach 
(right), the positron comes from t = +00 and can go forward and backward in time 
interacting a number of times with the ions. 
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5 Coulomb corrections for the heavy ion case, based on the Bethe- 
Maximon approach 

In the last section the exact solution of the one-particle Dirac equation in the 
high energy limit was discussed. Using the eikonal/ sudden approximation at 
high energies it was claimed in several publications that Coulomb corrections 
are not present. This claim must be wrong, since one has to reproduce the 
Bethe-Maximon results for photoproduction, see Sec. [2] in the limit where the 
charge of one nucleus Z\ is small, i.e. Z\<x <C 1. In Sec. |3]we dealt with pair 
production in the "usual" approach of QED, that is the perturbative calcula- 
tion and summation of Feynman diagrams. We saw that in the papers which 
we briefly reviewed very powerful techniques could be used to simplify the 
calculations at high energies. Also various terms could be summed analyti- 
cally. In the end an agreement with the Bethe-Maximon results was found. 
Yet, people not so familiar with these techniques might feel a little lost and 
other approaches (still for the high energy case) might be preferable. This is 
the subject of the present section. 

Before starting it may be mentioned that in Ref. [20J it was already realized 
that there was a problem with higher order Coulomb corrections, see Chap. 7.3 
there. There are corrections of the Bethe-Maximon type (class (ii) and (iii) 
in our notation) and corrections due to class (iv) (with n,n' > 1). An 'inter- 
polating' function f(Zi, Z 2 ) was conjectured in (20] (see Eq. (7.3.7)), in order 
to take care of these corrections but the question remained inconclusive. The 
convincing answer for the Coulomb corrections was first given in Ref. [36] . 



5.1 Coulomb corrections in the equivalent photon approximation: the ap- 
proach of Ivanov, Serbo and Schiller 



Let us now briefly review the approach of Ref. [36J . These authors start with a 
classification of the pair production amplitudes M nn i according to the number 
of photon lines attached to each ion, that is, the classification introduced in 
Sec. [TJ The sum of the amplitudes is written as J2nn' M nn > = M# orn + Mi + 
Mi + M 2 , corresponding to + + + M^) in our notation. The 
Born amplitude contains the one-photon exchange with each nucleus (class 
(i)), the amplitudes M 1 and Mi correspond to classes (ii) and (iii), M 2 are 
the processes of class (iv). In their Eq. (8) the total cross section is classified 
according to 

<? = O-Born + 0"l + + &2, (26) 



22 



where 



doBorn ~ \M Born \ 2 = |il%)| 2 (27) 



and 

dtn ~ 2Re(M Born M*) + |M!| 2 = 2Re{M {€) M* {ii) ) + |i% } | 2 (28) 



etc. Now the authors estimate the leading logarithms L = In (71 72) = ln(7 2 ) 
(see Eq. (1TB]) in Sec. [3]) appearing in the cross sections cxj. Due to the integration 
over the transferred momenta there are two large logarithms L in a Born ~ L 2 
and one large logarithm L for o\ ~ L and cfi ~ L. There is no large logarithm 
for cr 2 ~ L°. (Eventually, there will be another logarithm due to the integration 
over the energy of the e + e~pair.) Therefore this last (and most troublesome) 
term can be neglected, to a good accuracy, in the high energy limit when 
calculating the cross section. Please note that this argument is true for the 
cross section. As we will discuss below it does not hold for impact parameter 
dependent probabilities, especially for small impact parameters b. 

In order to study the higher order effects, the method of equivalent photons 
is used. Thus the higher order effects are reduced to a study of the higher 
order effects in the photoproduction (by real photons). These corrections are 
well known, see Sec. [2J and given in Eqs. (10)— (13) of Ref. [36J. Their final 
result is the expression Eqs. (24)-(27) for a\, and a corresponding one for 
6\. As discussed above, it is proportional to L 2 . Further, it is estimated that 
Gz/cTBorn ~ (Za) 2 /L 2 . Therefore this last class of processes will give a contri- 
bution so small that one can neglect this contribution for practical purposes. 
In order to see these effects either a very good accuracy would be needed, or 
one should find a way to increase their relative importance. One possibility 
to enhance the importance of small impact parameters is discussed in Sec. |HJ 
Their final result (Eq. (27)) for the Coulomb correction is then 

56 

o-coui ~ en + cr\ = — —a f(Z)L 2 , (29) 



where / is given in Eq. ffl3|) in Sec. [2] and o"o = - ^ 2 2 ■ 

It is also noted that the large transverse momenta p 2 ^> m 2 give a negligible 
contribution to the total cross section, but the experiments of STAR at RHIC 
(see Sec. [8] below) are just in this region, since one needs the high transverse 
momentum to detect the electrons. 
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5.2 How to regularize the exact solution of Baltz and McLerran and Segev 
and Wells: the approach of Lee and Milstein 



As shown above we know that there are indeed Coulomb corrections. They 
are predominantly given by the Coulomb corrections which are present in the 
photonuclear case. 

In order to see the Coulomb corrections in the exact solution of the Dirac equa- 
tion, we need to study in more detail the regularization that is needed in the 
sudden limit. In this limit all longitudinal momentum transfers are neglected. 
In this case the cross section would diverge, as large impact parameters con- 
tribute to arbitrary large values. Therefore a regularization of the result is 
needed. This is done in a certain way, by adding a longitudinal component 
to the final result in [29]. Another approach is to calculate only the deviation 
of the full results from the Born result, as it is only the Born result, which 
leads to a divergence. In [53,54] it is shown how to obtain reasonable results 
from the expression of Segev and Wells [55] and Baltz and McLerran |29j . The 
problem of the slowly decaying Coulomb potential is solved by introducing 
a screening of the Coulomb potential. Eventually the screening parameter is 
extended to infinity. The procedure of [51] and [53] is a nice example of the 
power of analytic methods. In order to outline their reasoning we would like 
to repeat some of the key steps of this work. For details we have to refer to 
their work. 

The authors of [52] start from an expression derived by Segev and Wells [55] 
and Baltz and McLerran [29] for the total pair production cross section, see 
their Eq. (1) 

* = t£^>i 2 ^® +a - WmW- (30) 



Here A; is a vector lying in the plane transverse to the momenta of the ion(s). 
The amplitude M is given by 



uto=m\ s - {k -^tT - 2 i- 

—p + q- — [k — p±) z — m z 
-a-(k- q ± ) + 7 m 

-j+]u(-q), (31) 



—p-q + — (k — q±) 2 — m? 



where p and e p (q and e q ) are momentum and energy of the electron (positron), 
u(p) and u{— q) are positive and negative energy Dirac spinors, a = 7°7, j± = 
7°±7^, 7 M are the Dirac matrices. We use light-cone variables p± = e p ±p z and 
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p+q-k 



q 



Fig. 7. One of the two lowest order diagrams for the pair production from two 
(virtual) photons with momenta k and p + q — k. 

q± = e q ± q z . The function F(A) describes the interaction with the Coulomb 
field of one of the ions. It is proportional to the electron eikonal scattering 
amplitude for one of the potentials V(r): 



For an unscreened Coulomb potential V(r) = —Za/r the integral for the 
Glauber phase x(p) becomes divergent and a regularization is required. One 
can choose, e.g., V = —Zaexp(—r/a)/r, i.e., the Coulomb potential is smoothly 
cut off at a radius a. Eventually, the radius a will tend to oo. It is shown in 
[53f54"] that the result is independent of the particular method used, as long 
as V{r) — * —Za/r for r — > 0. 

Using Eq. (9.6.23) of Ref. |56j and the substitution dz = rdr/y/r' 2 — p' 2 we 
find that x(p) = I dzV(b, z) = —2ZaK (p/a). Using the following well-known 
representation of the Bessel function (see e.g. [56J) 




(32) 



where 




(33) 



— oo 




(34) 



o 



one finds Eq. (11) of |54j : 




(35) 



o 
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Let us now take the limit a — > oo for a fixed A ^ 0. For p/a <C 1 we have 
Kq ~ In (p/a) and we get an integral of the type 



oo 

F(A) «27r hm / dppJ (pA)(exp2zZa:in(p/a) 



1). 



(36) 



This is a well-known generalization of the Bethe integral and one obtains 
T(l-iZa) ( 4 \ 

F = ™ Za WTIzJj{^) ■ < 37) 



Up to an overall constant phase, we see that the effect of the higher order 
Coulomb interaction consists of a change of the photon propagator 1 / A 2 to 
1/ (A 2 ) l ~ %Za ). This result was already found in |29j. Taking the square of F(A) 
one gets 



l^(A)| 



(A 



1 

A4> 



(38) 



independent of Za : which is the reason why the absence of Coulomb cor- 
rections was found in their case. But one should keep in mind, that A is the 
transverse part of the momentum transfer from the ions only In [53.54J a recipe 
is given that is used in order to obtain the "correct Coulomb correction". It 
is necessary to do the integration over k first in Eq. (130]) using the function 
F(k) with the regularized phase and then to remove the regularization only 
afterwards. 

Before the Coulomb correction to electron positron pair production itself is 
studied in refs. [531154] these authors studied the integral 

G = j^ 2 k\\F{k)?-\F\k)\% (39) 

where F°(A) = —i J d 2 pexp(—iA ■ p)x(p) corresponds to the Born expression. 
We integrate F(A) of Eq. (1351) by parts over p to obtain: 

oo 

Z7T1 r 

F(A) = —j dppJ 1 (Ap)x , (p)exp(-i X (p)). (40) 
o 

In order to obtain this result recursion relations for the Bessel functions are 
used (see Eq. (9.1.27) of |56j). The function F°(k) can be obtained from 
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Eq. (140]) by omitting the exponent in the integrand. Now the expression for G 
reads 



G = 2n J dkk J J dp 1 dp 2 piP2Ji(kp 1 )Ji(kp 2 )x(pi)x(P2) 
o oo 

x(exp(-zx(pi)+zx(p 2 ))-l). (41) 

A naive interchange of the order of integration would lead to G = 0, since due 
to the orthogonality relation for the Bessel functions the integration over k 
leads to a term proportional to S(pi — p 2 ). 

Instead the authors of Ref. [53] proceed in another way: they integrate k up 
to a finite fixed value Q, instead of oo. According to Eq. (11.3.30) of Ref. [5B] 
(for fM — v — 0) this gives their Eq. (7). 

The final result for G is then found to be 

G = -8tt (Z a) 2 f(Z a), (42) 



where, again, we encounter the ubiquitous function f(Za) of Eq. (ITHj) above. 

Now we proceed to the Coulomb correction itself. The Coulomb corrections 
are defined in the usual way discussed above as the difference of the full result 
to the Born result. One obtains three terms, see Eqs. (7) and (8) of [51]. In 
[53f54"] it is argued that the main contribution to the integrals comes from the 
region of small k. Accordingly, the matrix-element M of Eq. (|3T|) . is expanded 
as 

M = hk'jMij, (43) 



where k' = q + p — k. In [54] |Mjj| 2 is then calculated explicitly. This decom- 
position into three terms corresponds to the decomposition already discussed 
above in Eq. (126]) : 

O = (?Born + o\ + 0% + (J C A B = ° Born + 0\ + G\ + CT 2 - (44) 



We have a Bo m ~ L A L B , a c A ~ G A L B , cr B ~ L A G B , and o c AB ~ G A G B . The 
quantities G AtB are defined in Eq. (13"§|) and their value given in Eq. (T4"2~l) . The 
quantities L A B are defined by 



J A,B 



d 2 k 



k 2 \F%\ 2 . (45) 
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The necessary integrations are performed with logarithmic accuracy. The L 3 
behavior of the Born term is recovered. (Note that the full Racah result, see 
Eq. (1161) in Sec. [31 contains terms proportional to L 2 and L also). The Coulomb 
corrections a c A and a c B agree with those of [3B|, see Eq. (1211 . They show the 
typical L 2 = (ln7) 2 behavior. A new result is the last term, a AB , where the 
Coulomb correction applies to both nuclei. It scales only with L. This shows, 
as discussed above, that it is smaller by a factor of L 2 compared to the Born 
result and therefore is less important in the high energy limit. 

In a recent numerical calculation by Baltz et al. these results were checked 
by a calculation starting from the expression Eq. (1501) directly and doing the 
integrations numerically [5T|58II59II6U|I21| . We refer the reader to these papers 
for more details. 

Calculations using the exact solution as given in [55;29j were already done 
for small impact parameter. As discussed above the size relation between the 
different terms is only true for the cross section, but not for the impact param- 
eter dependent probability (see also Sec. E] below). It was found that Coulomb 
corrections are present at small impact parameter. This can reduce the prob- 
ability by up to 50%. These results were confirmed by a calculation in [6T]. 
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6 Higher Order Effects in Electron Pair Production: Multiple Pair 
Production and Bound-Free Pair Production 



In the previous sections we have mainly dealt with one special kind of higher 
order effects, the so-called Coulomb corrections. By definition they describe 
the effect of higher order Coulomb interactions on the one-pair production. 
They are well understood in terms of the Bethe-Maximon theory, or, equiva- 
lently, the approach of summing Feynman diagrams [27] . Higher order inter- 
actions can also lead to a second kind of processes: multiple pair production. 
Multiple pair production is the most important process of this type and will 
be discussed in this section. There will also be Coulomb corrections to mul- 
tiple pair production and we will see that they can be incorporated into the 
existing framework. 

In the high energy limit discussed in Sec. [5] the oppositely charged electrons 
and positrons enter in a symmetric way. This is in contrast to the second 
order correction for large angle pair production of [13], see the discussion in 
Sec. HJ An obvious asymmetry also occurs when the positively charged nucleus 
interacts with the electron to form a bound state (K-,L-,. . . shell capture). This 
will be dealt with in Sec. 16.21 We conclude with some remarks regarding muon 
pair production. In this case the finite size of the nucleus needs to be taken 
into account, which then leads to differences in the importance of Coulomb 
corrections and multiple pair production cross sections as compared to the 
e + e~ case. 



6. 1 Multiple pair production 

In heavy ion collisions with 77 = Zx ^f » 1 a classical impact parameter b 
can be defined. At the high ion energies considered here the classical path can 
be taken as a straight line, see e.g. Chap. 2 of [17] or [23] for details. One can 
then calculate impact parameter dependent probabilities for pair production 
P(b) in the semiclassical or Glauber approximations. The total cross section 
is given by the integration of this probability over all b. 

There is a handy formula for the impact parameter dependent (single) pair 
production probability in in lowest order, as shown in [27]. Quasireal 
photons from each nucleus collide to produce an e + e~-pair, the total cross 
section is denoted by <r 77 (s), where s is the square of the cm. energy. Using a 
simple version of the impact parameter dependent equivalent photon spectrum 
(valid for 1/m <C b <C 7/m) and the relation J^ a Y°>y( s ) — "if" ^2 one 
obtains 153,621 
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PW(6) 
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9tt 2 



m 
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28 ZfZ^a 



9tt 2 



2 In 7 z 2 a6 — 3 ln(m&) ln(m6) 

2 



for 1 < m& < ^ 



In 



for -fi ab < m& < 7 z 2 afe . 



(46) 



This equation displays nicely the dependence on the impact parameter b and 
the electron mass m. It also shows directly that for sufficiently large values of 
Zi, Z2 and 'jiab the first order probability P^(b) exceeds one, i.e., unitarity is 
violated. We further note that this calculation underestimates the probability 
in the range of small impact parameters b ~ A = — . The equivalent photon 
approximation used to derive this approximation is not justified in this region. 
More exact numerical calculations need to be done. 

The impact parameter dependent probability P^\b) in lowest order was cal- 
culated numerically in [63164] . The quantity P^ (b) found there is larger than 
the one given by Eq. ( 14*61) . and values larger than one are possible at RHIC 
and LHC. As an example P^(b = 0) = 3.9 (1.6), and pW(A) = 1.5 (0.6) [63] 
for LHC (RHIC) are found. 

The problem of the violation of the unitarity of P(b), that is, that P(b) will 
exceed unity for sufficiently large beam energies, was mentioned in [20]. It led 
to a series of studies starting around 1990. It was found that the production of 
multiple pairs in a single collision restores unitarity [SS]. e + e~-pair production 
at high energies can be studied in the sudden approximation, for not too large 
impact parameters and invariant masses. I.e., we can neglect the time-ordering 
operator in the Dyson-expansion of the S-matrix to get: 

S = exp(-i J V(t)dt). (47) 

—00 



One introduces the pair annihilation and creation operators bi and b\, where i 
denotes collectively the state of the pair (momenta and spin projection). To a 
good approximation the pairs can be treated as "quasibosons" with the usual 
commutation relations: 

[M] = <%. ( 48 ) 



We can then write 

00 

J V(t)dt = J2uibl+u*bi, (49) 



where Ui is a c-number which describes the probability amplitude for the 
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Fig. 8. Third order contribution to the one-pair production. 

production of a pair in state z in lowest order, see e.g. Sect. 7.1 of [20]. In 
principle, the interaction V also contains rescattering terms, i.e. terms, where 
the electron (or positron) of a pair changes its momentum and spin. We neglect 
these terms. We use the Baker- Campbell-Hausdorf identity 



,A+B 



e A e B e 



h[A,B] 



(50) 



which is valid for two operators A and B for which the commutator is a 
c-number. With this identity one can directly obtain an expression for the 
production of multiple pairs with given quantum numbers, valid to all orders. 
In Fig. El we show a third order contribution to the one-pair production. The 
other e + e~ pair which is created is annihilated (by a "light-by-light scattering 
diagram"). This is quite similar to the emission of multiple soft photons due 
to a classical current, see e.g. [12]. One can say that a coherent state [66] of 
e + e~ pairs is produced. 

Similar approaches were followed in [67] and [65]. Another early work in this 
direction was given in [6^170] . All studies essentially make use of the "quasibo- 
son" approximation, see Eq. (1501 above. It is then found that the probability 
to produce N pairs P(N, b) is given by a Poisson distribution: 

P(N, b) = exp [-PW (&)] , (51) 



where P^\b), the probability for pair production calculated in lowest order 
(which can become larger than one). The quantity P^'(b) acquires the mean- 
ing of the "average number of e + e~ pairs" produced in a single ion collision: 

(N(b)) = NP{N, b) = PW(6). (52) 

JV 



Another approach without the use of the quasiboson approximation was done 
in [51] and later also in [TTJ. There it was found that quite generally the 
amplitude for iV pair production by an external field is given by 

S N =(0\S\0)^sgn(a)stj, 1 ---st; l , N) , (53) 
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Fig. 9. Graphical illustration of the TV-pair production process. The interaction 
with the external field is shown as crosses. The production of a pair is described 
by a fermion line coming from and leaving to the future, interacting an arbitrary 
number of times with the external field. The vacuum- vacuum amplitude (0| S\0) 
corresponds to the sum of all closed fermion loops. 



where ki,k are the quantum numbers (momenta and spin projection) of elec- 
tron and positron, respectively, and a denotes a permutation of {1, • • • , N}, 
see also Fig. [91 The vacuum-to-vacuum amplitude is present in all QED cal- 
culations, but is often dropped if it is of absolute value 1. Here this amplitude 
is < 1 and is due exactly to the contribution of all light-to-light diagrams, 
see Fig. |HJ The amplitude s + ~ is the one for the production of a single pair, 
corresponding to a single fermion line. Neglecting the antisymmetrisation of 
the amplitudes, one gets the N pair amplitude as the product of the vacuum- 
vacuum amplitude and the single-pair production amplitudes. From this one 
gets again a Poisson distribution for the iV-pair production probability. As 
before the approximation introduced lead to an independence of each pair 
production process from all others, which is the essential ingredient needed to 
get the Poisson distribution. For further details we refer to Chap. 6 of [TS] . 

With the impact parameter dependent probability P^'(b) calculated, e.g., in 
[63|M] one can use Eq. floTl) to obtain the probabilities for iV-pair production 
P(N,b), see Fig. [TUJ One can see that for impact parameters b ~ 2R up to 
about X on the average 3-4 pairs will be produced in Pb-Pb collisions at 
the LHC This means that each photon-photon event — especially those with 
high invariant mass, which occur predominantly at impact parameters close to 
b > 2R — is accompanied by the production of several (low-energy) e + e~ pairs 
(most of them however will remain unobserved experimentally). Integrating 
over the impact parameter the total multiple pair production cross section 
was given in [72] , see Fig. [TTJ 
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Fig. 10. The impact parameter dependent probability to produce ./V e + e~-pairs 
(N = 1,2,3,4) in one collision is shown for RHIC (7 = 100, Au-Au) (a) and for 
the LHC (7 = 2950,Pb-Pb) (b). Also shown is the total probability to produce at 
least one e + e~-pair J2n=i b) = 1 — P(0, 6). One sees that at small impact 
parameters multiple pair production dominates over single pair production. 

10 s 

H) 2 

^ 10' 
X 10° 



Iff 2 
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Y 

Fig. 11. Total cross section for the production of one, two and three pairs in PbPb 
collisions is shown as a function of the ion energy. 

The Poisson distribution is also valid if, e.g., Coulomb corrections are included 
in the calculation of the single pair production. As seen in Eq. (153]) they need to 
be taken into account for a single lepton line to get a -P 1 (6) including Coulomb 
corrections. This -P 1 (6) was calculated for the high energy limit, as discussed 
in Sec. EJ in [73l6Uf61j . It was found that the Coulomb correction reduce -P 1 (6) 
up to 50%. The single pair production probability P(b) falls off essentially as 
1/b 2 . The total cross section is therefore dominated by large impact parameters 
and not very sensitive to the multiple pair production effects at small b. In 
[62] the reduction of the exclusive one-pair production cross section a(N = 1) 
was estimated to be -6.4% for RHIC and -4.7% for LHC, using the numerical 
results of [51172] . 









— —""* N=l 








W 












Multiple pair production cross sections are dominated by impact parameters 
around b m A. It would be of interest to establish this effect experimentally. 
In addition, Coulomb corrections are enhanced for these close collisions. At 
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Fig. 12. Bound free pair production is the pair production process, where the elec- 
tron is produced into a bound state of one of the ions. 

the SPS (CERN) the effect was looked for, see [71], but only an upper bound 
could be given, which is still above the theoretical prediction. The possibilities 
at the LHC will be briefly discussed in |HJ 

6.2 Bound- free pair production 

Bound-free pair production is the pair production process, where the electron 
is created in a bound state with one of the ions, say Z\ (see Fig. [TJj) 

Zl + Z 2 -> {Z x + e~) KjL> ... + e+ + Z 2 . (54) 

Since our previous review [18J essentially no further calculations of bound- 
free pair production were made. The most careful state-of-the-art calculations 
were done in [12]. It is a Born-approximation calculation, using Dirac wave 
functions for the electrons and positrons in the field of a point nucleus with 
charge Z\. In this approximation the cross-section scales with Z\. At the high 
energies relevant for RHIC and LHC this is a good approximation, see also 
Sec. Higher order processes of the class (iv) would lead to a deviation from 
the Z\ scaling. They have been calculated using the ultrarelativistic limit of 
the Dirac equation in [75J. Only a small reduction of the order of 1% was 
found. 

In an "Einstein special issue" electron-positron pair production in relativistic 
heavy ion-atom collisions was reviewed in [76J in a general way, not restricted 
to the high energy limit. In this paper a discussion of the bound-free pair pro- 
duction mechanisms is also given. While at the high energies the "excitation 
type" mechanism studied in [12] and in other papers dominates, the "trans- 
fer type" mechanism is also important at lower beam energies. "Excitation" 
here means, that the electron and the positron are best described by wave 
functions in the target Coulomb field, "transfer" means that the positron is 
in the Coulomb field of the target, whereas the electron is "transferred" to 
the Coulomb field of the projectile. For the lower beam energies it becomes 



34 



more proper to treat pair production by a two-center calculation. At high en- 
ergies, 'class (iv)' is relatively unimportant, see above. On the other hand, at 
lower energies this type of process becomes increasingly important and causes 
further calculational complexities. 

The cross-section for bound-free pair production at high energies rises only 
with In 7, i.e., 

a = Ah\^ + B. (55) 



It is only a small fraction of the total cross section of free pair production, 
which rises with (ln7) 3 . The parameters A and B have been given for various 
systems in [12]. Nevertheless, this process is an important one for the operation 
of the Pb-Pb beam at LHC: the ions which have captured an electron from 
the vacuum have a charge-to-mass-ratio different from the beam particles and 
they will hit the wall at a rather well defined 'hot spot'. This is discussed 
further in Sec. [8] below. 



6. 3 Higher order corrections to muon pair production 



In this review we focus on the production of electron-positron pairs. Rather 
similar phenomena will appear in muon pair production, scaled with the corre- 
sponding mass. The larger muon mass has an effect on the different processes: 
in [27] it is remarked that the nucleus is an extended object with the inverse 
radius 

A = -1 ~ SOMeV. (56) 
R 

At such scales the electric field differs substantially from the pure Coulomb 
behavior. This is especially of interest for muon pair production, where we 
have the inequality 

m < A < m M . (57) 



This means that a simple scaling of the cross section by m 2 /m 2 ^ is not valid. 
It is shown in [27J that there is a suppression factor of A/m^ for additional 
photon exchanges. This is because the form factor cuts out the short wave- 
length photons with k > A and the remaining long- wavelength photons cannot 
resolve the dipole formed by the lepton pair which is of extension ~ 
This leads to a reduction of the Coulomb effects for the muon case compared 
to the electron case. On the other hand unitarity corrections, due to the ad- 
ditional production of electron-positron pairs together with the muon-pair, is 
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enhanced. This can be understood as the impact parameter range in the muon 
pair case is much smaller than in the electron case and the probability to pro- 
duce (additional) electron-positron pairs is then large. For further discussion 
we refer to [27]. 
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7 Transition from the adiabatic to the sudden regime 



The physics of e + e pair production in fast (v ~ c) and slow (v ~ 0) heavy ion 
collisions (Z unit cd atom ^ 137(173)) is very different: for the slow collisions there 
is spontaneous pair production in supercritical fields, a non-perturbative effect. 
For fast collisions we can apply perturbation theory. Two-photon production 
dominates, as was discussed above. Let us first mention the corresponding 
situation of atomic ionization and electron-positron-pair production in a time- 
varying spatially-constant electric field of the form 

E(t) = Eocosut. (58) 



The problem of ionization of atoms in a time varying electric field is very 
similar to the problem of e + e~ pair production, which may be viewed as 
ionization of the negative energy Dirac sea, see Sec. HI For ionization there is 
a similar transition from a slowly varying field to the high frequency case. The 
ionization of an H-atom in a time varying electric field is a textbook problem, 
see, e.g., p. 739ff, problem XVII in [77J. One limit — large frequency u and a 
small field strength E — is the usual photo-effect. The other limit is ionization 
in a static electric field. For a sufficiently strong static (u = 0) electric field, 
of the order of E a = ^^r- = E c a 3 , the electron disappears classically into the 
continuum; for lower field strength it can tunnel through the barrier, see, e.g., 
[78]. 

The critical atomic field strength E a is given by the condition eE a dB = Ry, i.e., 
the work done in the field over the typical length scale (the Bohr radius ag) is 
then equal to the typical energy scale (the ionization energy Ry = e 2 /(2a#)). 
(To obtain the critical field strength for the Schwinger case, discussed in Sec. HJ 
one needs to replace the Bohr radius by the Compton wavelength of the elec- 
tron A = h/(mc) and the Rydberg energy by the electron rest mass mc 2 .) The 
field strength occurring in current lasers, see, e.g., [79] is strong enough to lead 
to such ionization processes: in Chap. VI of [79J the optical field ionization 
of atoms is reviewed. The important parameter which divides the strong field 
and the weak field limits is the Keldysh parameter 7^ (see Eq. (21) of [79J). 
For 

1/7* - ^ « 1 (59) 



we have the perturbative limit, for 

1/7* > 1 (60) 
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we have the strong field regime. We do not go further into this but turn to 
e + e~ pair production case, which was studied in [80J. 

The adiabaticity parameter (see also [T3p^] ) is given by Eq. (47) of [80J: 

mcu} fai\ 

IB! = (61) 



It corresponds to the Keldysh parameter jk of the atomic physics case. The 
Compton wave length of the electron X, now replaces the Bohr radius ag. Their 
final formula for the pair-production probability per unit time and volume is 

aE 2 ixm 2 

w = o~T7 \ I 1 /n 77 vV ex P ^g{lBi))- 62 

2ir(g(iBi) + 1/2^' (7b/)) eE 



The function g(jsi) is given by Eq. (45) of [15], it interpolates between the 
two extreme situations. In the limit jbi *C 1 it is given by 

9(lBi) = l-llBi + 0(j BI ) (63) 



and for jbi 1 they find 

g(lBi) = — ln(2 7s/ ) + 0(1/ 7 b/). (64) 



The parameter jbi describes the transition from the high-field low- frequency 
limit (tbj <C 1) to the low-field perturbative regime Tbj S> 1). In the jbi ^ 1 
limit, one obtains the Schwinger formula, see Sec. [TJ For j B i ^ 1 the pertur- 
bative result (see also Eq. (11) of [MUl]) is found: 

9 



d n c fe E\ 



d 3 xdt 



4^UJ • (65) 



It corresponds to n-th order perturbation theory, where n is the minimum 

2mc 2 
Ttuj 



number of quanta required to create an e + e pair (n > ^jp 2 ^ ^> 1). 



For the heavy ion scattering the space-time dependence of the electric fields 
is much more complicated and one cannot obtain such instructive analytic 
results. Let us roughly estimate a parameter jhi for heavy ion collision corre- 
sponding to ■ysi- we put oj max ~ 1/At = With the maximum field strength 
E ~ -gp-, see Eq. (j3J) in Sec. [fl we obtain: 

mcbv p e i 

lHI = ^ = Ap' (66) 



38 



where p e i = mc is a typical electron momentum scale and Ap = corre- 
sponds to the momentum transfer in a Coulomb collision(see Eq. (jlj) in Sec. [I]). 
For a typical impact parameter b one can take the Compton wave length of 
the electron, i.e., b rs H/(mc) . So we have 

lM = v/(Za). (67) 



In the relativistic case, v ~ c, this is > 1, i.e., we are in the low-field high- 
frequency limit. For v <C 1 we are in the strong-field non-perturbative limit. 

Our conclusion from this chapter is the following: The slow collisions, with 
their static overcritical fields are in a different regime of the adiabaticity pa- 
rameter as compared to the fast (ultra-) relativistic heavy ion collisions: the 
perturbative approach, (sometimes up to spectacularly high orders, see [23]) is 
well justified. The crossover between these two regimes can then be expected 
to occur in the region where the electron velocity is of the order of the velocity 
of the ions, i.e. v ~ Za, which corresponds to 7 « 1 ... 2. 

There is strong support for a future x-ray-free-electron laser at Desy (Ham- 
burg), see [81] and one can expect exciting new results. The possibilities of 
optics in the relativistic regime, which are opened up by ultraintense laser 
pulses, are described in [82j. We mention the talk by A. Ringwald [83] or [13] 
and the theoretical calculations of pair creation in the fields of an X-ray free 
electron laser [51]. On the other hand, from what we said above, we have per- 
turbative physics from the heavy ion collisions: the fields get very strong, with 
increasing ion energy, but, at the same time, the interaction time decreases in 
the same way. 
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8 Comparison to experiment and outlook to LHC 

Since our last review of this subject in [18] the heavy ion collider RHIC has 
come into operation. Both the STAR and the PHENIX collaborations have 
published measurements on e + e~ pair production. In the case of STAR the 
data was taken for events accompanied by nuclear breakup. In this way one 
is able to measure pair production at small impact parameters. Due to the 
high Lorentz factor 7, the electromagnetic fields are stronger than at SPS 
or AGS energies. For a discussion of experiments at these accelerators we 
refer to previous reviews |18I|85II17] . As explained in Sec. [7] we are now in 
the perturbative limit. Also multiple pairs and bound free pair production 
processes can be measured in principle. It is the purpose of this section to give 
a short overview of these results and to compare them with the theoretical 
analysis, especially to the one done in [SB]. Bound- free pair production was 
studied at RHIC as well. In the end we give a short outlook to the possibilities 
at the LHC. 



8.1 e + e pair production in ultraperipheral collisions at RHIC 

The STAR collaboration has published data on e + e~ pair production accom- 
panied by nuclear breakup in ultraperipheral Au-Au collisions at a center of 
mass energy of 200 GeV per nucleon [35] • This process is illustrated in Fig. [T31 
The nuclear breakup of the gold ions is predominantly due to the electromag- 
netic excitation of the giant dipole resonance (GDR), see [X8"fTT] . The breakup 
neutrons are detected in the ZDC (Zero Degree Calorimeter). It can be shown, 
that the amplitudes for these processes factorize [23J and that the cross sec- 
tion can be calculated by integrating the product of the probabilities over all 
impact parameters 



where PcDRib) denotes the probability of exciting one ion and (Pgdr(P)) 2 
is the probability to excite both ions. In this way pair production at small 
impact parameter can be selected. The average impact parameter is given 
approximately by [23] 




(68) 



b 



Jd 2 bbP(b) 8R a 



19 fm, 



(69) 



Jd 2 bP(b) 3 
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Fig. 13. At STAR pair production in connection with the electromagnetic excitation 
of one or both ions was measured. 

The impact parameter dependent probability was calculated in lowest order in 
[8"6] based on the method developed in [53|87f72j . Only electrons and positrons 
with transverse momenta of p t > 60MeV/c with rapidity \y\ < 1.15 could be 
detected. This limited detector acceptance is taken into account in the calcula- 
tion. Within the experimental accuracy a good agreement was found between 
the experimental and theoretical results not only for the total cross section 
but also for the distribution in energy and transverse momentum of electron 
and positron, see Figs. HU [15] and [16] and in reference [381188] . The results 
of the calculation agree also with a simpler calculation using the equivalent 
photon approximation in most cases. The exception is the distribution of the 
transverse momentum of the pair. 

In the lowest order calculation electrons and positrons have an identical dif- 
ferential distribution. No asymmetry was also seen in the experimental data, 
which has a large statistical error. Also the total cross section was found to 
be in agreement with the lowest order calculation. No higher order reduction, 
as predicted by the Bethe-Maximon theory was observed, which should be a 
reduction of about 17%. But this is mainly due to the restriction on the phase 
space by the experiment, where these corrections are not expected. 



8.2 Bound-free pair production at RHIC and the forthcoming LHC(Pb-Pb) 

The bound-free pair production (BFPP) (see the discussion in Sec. 16.21) is 
a process which is of practical importance for the operation of the LHC in 
the Pb-Pb mode. An electron-positron pair is produced, where the electron is 
produced not as a free particle but into a bound atomic state of one of the 
ions, 

Zt + Z 2 -> (Z 1 + e-) 1Sl/2i ... + e + + Z 2 . (70) 
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STAR eV Differential Cross-section vs. Pair Invariant Mass 
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Fig. 14. Comparison of the distribution of the pair mass for the STAR experiment 
to the theoretical QED calculations of [86] and to a calculation using the equivalent 
photon approximation (EPA). Taken from |38j 



STAR eV Differential Cross-section vs. Pair Transverse Momentum 
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Fig. 15. Comparison of the transverse momentum distribution of the electron and 
the positron for the STAR experiment to the QED and the EPA calculation. Same 
as in Fig. [T4l Taken from [38] 

As this changes the charge state of the nucleus, it is lost from the beam in 
the collider. Together with the electromagnetic dissociation of the nuclei (see 
[T8][rT] these two processes are the dominant loss processes for heavy ion col- 
liders. It has been realized [89], (see also [90|9T] ). that the BFPP process can 
also result in a localized beam-pipe heating: the atomic states are produced 
with a small perpendicular momentum of the order of m, and therefore leads 
to a narrow singly-charged ion beam (The process of electromagnetic disso- 
ciation is less severe in this respect as the momentum transfer to the ions is 
more spread out and the ions are not so strongly focussed on a single spot). 
These beams with altered magnetic rigidity will deposit their energy in a lo- 
calized region of the beam pipe and cause a localized heating. This can lead to 
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STAR eV Differential Cross-section vs. Pair Transverse Momentum 
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Fig. 16. Comparison of the transverse momentum distribution of the pair for the 
STAR experiment to the QED and the EPA calculation. Same as in Fig. [Fil Taken 
from |38j. 



the quenching of the superconducting magnets. The energy deposited per unit 
time is of the order of Watts. This limits the luminosity of the Pb-Pb collider 
at LHC [90J. The effect of BFPP is studied more quantitatively through sim- 
ulations of the ion showers in a dipole magnet in [92.93J. These calculations 
are based on the cross section given in [42] . In this reference, a comparison is 
also made with the results of other calculations. The authors of [92] conclude 
that ". . . the secondary beam of ions emerging from each collision point at 
the nominal peak luminosity L = 1.0 x 10 27 cm _2 s _1 is not likely to quench a 
dipole magnet". 

Due to the lower beam energy at RHIC, the energy deposit there is smaller and 
therefore not important as a limiting factor to the beam luminosity. However, 
it can serve as a testing ground for the theoretical numbers given in [12] ■ This 
is especially important, as the experimental result with the highest energies is 
from a fixed target measurement with an energy of 160 GeV per nucleon |94j . 
corresponding to a 7 = 9.3, still far away from the situation at the LHC. 

First results of a measurement of ion beam losses due to BFPP in Cu-Cu 
collisions at RHIC were reported in [95]. Due to the small number of events 
observed the authors were not able to give a value for the cross section, but 
they conclude that within the uncertainties associated with the experimental 
conditions the data are consistent with the signal expected from BFPP. It is 
a valuable test of the ability to quantitatively predict the BFPP effect for the 
LHC. 
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8.3 Outlook on LHC 



In [86] also an outlook to the situation for the LHC is given. One of the 
drawbacks of all possible measurements at high energies is the fact that pair 
production is predominantly produced in the forward direction, whereas the 
detectors are designed towards central events. But most of the interesting 
effects, like Coulomb correction and multiple pair production will occur at 
small transverse momenta. In [96J it is discussed how the very inner part of 
the detectors can be used for pair production detection. The limits of the 
transverse momentum can be reduced to p t > 2.6MeV, which is comparable 
to the rest mass. This is explored in more detail in [SB]. It is found that about 
10% of all pairs produced together with nuclear excitation will come from 
multiple pair production. 

Another approach has been studied in [97]. Here the very forward CASTOR 
detector is used to measure pair production. Calculations in lowest order are 
done and it is proposed to use this as a luminosity measurement for the LHC. 
Of course the accuracy of such a luminosity measurement relies on the precise 
theoretical knowledge of this process especially in this region and how large 
Coulomb correction contribute to it. In combination with the nuclear excita- 
tion of the ions it would also allow to study pair production at small impact 
parameters. 
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9 Conclusion 



In April 1990 a workshop took place in Brookhaven with the title 'Can RHIC 
be used to test QED?' [98]. We think that after about 17 years the answer to 
this question is 'no'. However, many theorists were motivated to deal with this 
topic. The gradual progress which was sometimes quite tortuous, is described 
in this report. 

In this review we studied electron-positron pair production in heavy ion col- 
lisions at relativistic energies. There were quite a few papers in the last years 
approaching this topic from various aspects. Whereas the lowest order Born 
result was known for a long time, higher order (Coulomb) corrections have 
been under debate for the last decades. These corrections were studied using 
various methods. One of these methods follow the approach of Bethe and Max- 
imon, which was used to study Coulomb corrections to the pair production 
in photon- nucleus interactions (the Bethe-Heitler process). QED perturbation 
theory was also used and it was shown [27] that they yield the same result 
as the approach of Bethe and Maximon. Difficulties with obtaining analytical 
solutions appear when more than one photon are attached to each nucleus, 
and no firm results have been obtained up to now. However, for the total cross 
section these effects are small. 

It has been found in theoretical calculations that higher order effects can be 
quite sizeable for collisions where small impact parameters are selected. But 
an experimental confirmation has not been made up to now. 

The fields occurring in ultrarelativistic heavy ion collisions are very strong. 
So one may expect that some new kind of nonperturbative effects could arise, 
similar to the Schwinger mechanism for static electric fields. But the fields 
in the relativistic heavy ion collision act only for a short time and higher 
order perturbation theory is appropriate. We also discuss the problem of the 
transition from slow to fast collisions. However, this chapter is far from being 
solved. A new higher order effect in relativistic heavy ion reactions is the 
emission of multiple pairs in a single collision. This is quite a spectacular effect 
and it would be of some interest to observe it experimentally in forthcoming 
experiments at the LHC 

In addition to the intrinsic interest of the theoretical study of pair production 
at high energies there is the practical importance of these processes in the more 
general study of ultraperipheral processes. Pair production is a background 
process, e.g., to vector meson photoproduction. The very large cross section 
for bound-free pair production limits the maximum Pb-Pb luminosity at the 
forthcoming LHC. Thus it is mandatory to understand these QED effects to 
a high degree of accuracy. 
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Finally, we review the experimental results on pair production at RHIC and 
their theoretical interpretation. Results on the bound-free pair production 
mechanism obtained recently at RHIC are also discussed. Agreement with 
theory is good, and we expect a similar good agreement at the forthcoming 
LHC. This good agreement is of course also the reason why electron and muon 
pairs produced in ultraperipheral collisions may well be used to monitor the 
luminosity at colliders. This should not be surprising since QED is the best- 
tested theory we have. 
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